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ABSTRACT 

A sample of 59 sight lines to reddened Galactic OB stars was examined for correlations of the strength 
of the CO Fourth Positive [A^Ii — A^I]+) absorption band system with the ultraviolet interstellar 
extinction curve parameters. We used archival high-dispersion NEWSIPS lUE spectra to measure the 
CO absorption for comparison to parametric fits of the extinction curves from the literature. A strong 
correlation with the non-linear far-UV curvature term was found with greater absorption, normalized to 
E{B — V), being associated with more curvature. A weaker trend with the linear extinction term was 
also found. Mechanisms for enhancing CO in dust environments exhibiting high non-linear curvature 
are discussed. 

Subject headings: Dust, extinction — ISM: abundances — ISM: molecules — ultraviolet: ISM 



I. INTRODUCTION 

The extinction of starlight in the ultraviolet (UV; A < 
. 3200 A) can be characterized by three main components: 
a linear rise, a Lorentzian-like bump centered on A ~ 2175 
A, and a far-ultraviolet (A < 1700 A) non-linear rise. The 
general features of the extinction curve can be understood 
in terms of dust grain populations (Mathis, Rumpl, & 
Nordsieck 1977; Draine & Lee 1984, e.g.). Cardelh, Clay- 
ton, & Mathis (1988, 1989) showed that the overall shape 
, of the UV extinction curve can be estimated by the single 
' parameter Rv[= A{V) / E{B — V)]. However, individual 
' sight lines exhibit large deviations from the average ex- 
, tinction curve, particularly in the far-UV rise (Mathis & 

■ Cardelli 1992). A detailed understanding of these devia- 
\ tions can provide information about the nature of the par- 
. tides responsible for the extinction, putting constraints on 
' interstellar grain models. 

A number of studies have investigated the possible cor- 

■ relations between the various components of the UV ex- 
' tinction curve and gas phase abundances (Joseph et al. 

1989; Jenniskens et al. 1992), infrared emissions (Cox & 
Leene 1987; Hackwell et al. 1991; Boulanger, Prevot, & 
Cry 1994), very broadband structure (Jenniskens 1994), 
and diffuse interstellar band (DIB) absorptions (Desert, 
Jenniskens, & Dennefeld 1995). The Jenniskens et al. 
(1992) study was the first to positively identify a corre- 
lation of the strength of the non-linear rise with the abun- 
dance of a gas phase atom or molecule, finding that a 
larger column of CH implies a stronger far-UV non-linear 
rise, and the inverse relationship for CH'^. 

No such study has to date been performed for CO, a 
molecule of great astrophysical importance. Being the sec- 
ond most abundant molecule after H2 , it is used as a tracer 
of H2 gas and an indicator of the total mass of molec- 
ular clouds (Scoville et al. 1987; Bloemen 1989; Hunter 
et al. 1994). Having allowed rotational transitions that 
are readily observed in emission at radio wavelengths, the 
CO molecule is commonly used to probe dark molecular 
clouds. CO also has an extensive band system in the ul- 
traviolet, which can be used to investigate the molecular 



content of diffuse and translucent clouds when observed in 
absorption toward background stars. 

The wavelength coverage of the short-wavelength prime 
(SWP) camera of the International Ultraviolet Explorer 
(lUE) satellite, 1150 - 2000 A, spans as many as 19 bands 
of the Fourth Positive (A^H - X'^Yj+) [v'-Q) band system 
beginning with the (0-0) band at 1544.5 A and continuing 
to shorter wavelengths. For this study, we searched the 
lUE archival data for ultraviolet absorption of CO along 
59 lines of sight toward Galactic O and B stars for which 
UV extinction curves have already been published, and pa- 
rameterized as described below, to investigate the possible 
correlations with the various curve features. 

Fitzpatrick & Massa (1990, FM hereafter) have pro- 
posed a parameterization that accurately fits the extinc- 
tion curves with a single analytical expression using six 
parameters. Using this method, the extinction curve, nor- 
malized to unit E{B — V) and with x = A~^, can be ex- 
pressed as 

^(^) ^ — T7T = ci+C2X + czD{x;^,XQ) + CiF{x) (1) 



E{B - V) 



where 



D{x;-f,xo) 



(2) 



F{x) 



is the Lorentzian-like "Drude" profile, representing the 
2175 A bump with xq the bump peak and 7 its width, 
and 

0.53(a; - 5.9)2 -f 0.05446(a;- 5.9)3 ^ > b.^^im- 
X < 5.9/im" 

(3) 

is the far-UV curvature term, whose strength is determined 
by a single parameter, C4. 

2. THE SAMPLE 

Our sample consists of 59 reddened stars (0.24 < E{B — 
V) < 1.09). The UV extinction curves for these stars have 
been previously determined using the pair method by FM 
and Aiello et al. (1988). The extinction curve parameters 
for stars from the Aiello et al. sample are taken from Jen- 
niskens & Greenberg (1993, JG hereafter), who used the 
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FM fitting routine. For stars that appear in both sam- 
ples, we adopt the FM numbers. The parameters from JG 
show a small (1-(t) systematic offset from the FM data, 
producing slightly smaller ci values and larger C4 values, 
approximately —0.1 and +0.1 respectively. These system- 
atic errors do not significantly affect the results of the anal- 
ysis. 

Of the stars in the FM and JG samples, only those for 
which high-dispersion short-wavelength lUE spectra exist 
are considered here. We also consider only stars of spectral 
type B5 and earlier and luminosity classes III and V, since 
spectral mismatch using the pair method can strongly af- 
fect the extinction curve fitting at later types and higher 
luminosities (Massa et al. 1983). Although these stars 
do not represent a purely random sampling of interstel- 
lar sight lines, they do sample different interstellar envi- 
ronments, including dense and diffuse media, as well as 
regions of recent early-type star formation. 

3. ANALYSIS 

For each star high-dispersion short-wavelength New 
Spectral Image Processing System (NEWSIPS) lUE spec- 
tra were obtained from the archive maintained by the As- 
trophysics Data Facility (ADF). The quality of the data 
varies strongly from object to object due to the differences 
in exposure time, brightness of the source, and number 
of spectra available. To improve the signal-to-noise per 
source, multiple spectra were co-added by weighted mean 
when available. However, no further processing was per- 
formed on the data. Figure 1 shows examples of both 
strong and weak CO absorption. 

We measured the equivalent width of the (2-0) band in 
each spectrum. A linear continuum was fit and the equiv- 
alent width measured by integrating over the absorption 
line profile (typically a span of 0.4 - 1 A). The i-cr error 
was computed by summing the normalized errors for each 
pixel in quadrature. For non-detections, the upper limit 
was taken to be the error computed across 0.4 A. The mea- 
sured equivalent widths and the lUE exposure numbers for 
each star are listed in Table 1. 

The (2-0) band, located near 1477.6 A, was chosen for 
several reasons. It is well separated from the nearby ^'^CO 
band at 1478.8 A and is relatively uncontaminated by 
other interstellar or stellar features for all the spectral 
types considered in this study. It is also the least sus- 
ceptible to background subtraction errors known to occur 
in NEWSIPS data because it lies near the sensitivity peak 
of the echelle spectrograph (Massa et al. 1998). It also has 
the highest band oscillator strength of the Fourth Positive 
band system and thus is the most likely to be observed in 
low SjN data. 

The resolution of the lUE echelle spectrograph (i? k, 
10'') combined with the Hmiting SjN (« 20 : 1) of the 
SEC vidicon detectors (Cowie & Songaila 1986) conspires 
to restrict our typical 2-ct equivalent width detection limit 
to ~ 30 mA. As wc will show in the Appendix, even at 
this equivalent width, all absorption lines are affected by 
saturation and thus the conversion of the molecular band 
equivalent width to a corresponding column density de- 
pends on assumptions of the gas rotational excitation tem- 
perature, T, and doppler velocity, 6. Although the various 
sight lines are likely to have different values for T and 6, on 
average the values for these parameters lie around T « 4 



K and 6 « 1 km s^^ (Gredcl, van Dishoeck, & Black 1991), 
as observed in higher quality data. We will adopt these 
values to convert the equivalent width to column density 
for the comparison to previous work and allow a statistical 
assessment of the robustness of our result with respect to 
the phenomenologically more meaningful column density. 

For sight lines whose CO absorption has been previously 
studied our measurements agree very well with the pub- 
lished data. In particular, the equivalent widths for the 
(2-0) bands of HD 21483, HD 47129, and HD 149757 of 
322 ± 27 (Joseph et al. 1986), 46 ± 17 (Tarafdar & Kr- 
ishna Swamy 1982), and 74 ± 7 (Wannier et al. 1982) all 
agree within the error with our measured values. Four 
more sight lines (HD 34078, HD 37903, HD 147933, and 
HD 193322) have been observed in other bands, either 
higher A — X bands or the shorter wavelength B — X and 
C — X bands. Column densities for all of these stars are 
compiled in Federman et al. (1994). By converting our 
measured equivalent width to column density as described 
in the Appendix we directly compare our results for all the 
aforementioned sight lines in Figure 2. 

For the purpose of comparison to the extinction curve 
parameters as specified in Equation 1, the equivalent width 
is normalized to unit reddening. This produces a measure 
of the relative abundance of CO to the total amount of 
matter along the line of sight. Figure 3 shows the correla- 
tion plots of this quantity versus each of the six extinction 
curve parameters for the stars listed in Table 1. 

These data were analyzed using the Astronomy SUR- 
Vival Analysis (ASURV) software package Rev 1.2 (Isobe 
& Feigelson 1990; LaValley, Isobe, & Feigelson 1992), 
which implements the survival analysis methods presented 
in Feigelson & Nelson (1985) and incorporates the up- 
per limits in a statistically appropriate manner. For each 
parameter of the extinction curve, the Spearman's Rho 
value, a non-parametric rank-order correlation coefficient, 
was computed. This value, which can range between 1 
and —1 for a perfect correlation and anti-correlation re- 
spectively, was tested for its significance by calculating a 
t-value, which is distributed approximately as Student's 
distribution with N — 2 degrees of freedom and is effec- 
tively a measure of the standard deviations from the null 
hypothesis. This allows the probability that the correla- 
tion could be drawn randomly from this distribution to be 
determined. These results are tabulated in Table 2. 

4. RESULTS AND DISCUSSION 

No significant correlations were seen with the extinction 
curve parameters xq, 7, or C3 (see Figure 3). The relation- 
ship with C2 (middle right panel of Figure 3), while not 
showing a linear correlation does show an interesting bi- 
furcation. (FM note that ci and C2 are strongly correlated 
and thus only the relationship with C2 is considered here.) 
Although a high C2 does not guarantee a large CO abun- 
dance, the largest abundances are seen when C2 > 0.5. 
Of notable exception is HD 147889. This highly reddened 
star {E{B - V) = 1.09) has a C2 of 0.151 but the ab- 
sorption suggests a column density as high as 3 x lO"'^ 
cm^^. However, it must be noted that this star has a 
high C4 value of 0.709. As C2 is inversely proportional to 
Rv (Fitzpatrick 1999), our result would suggest that lower 
CO abundances would be observed along sight lines with 
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Rv}!L 3-6. Cardelli (1988) noted a decline in the abun- 
dances of H2 and CH with increasing Ry and a similar 
break in the CH abundance at Rv= 3.5 for clouds with 
Av < 2. The interstellar average Rv is 3.1, which corre- 
sponds to a C2 of 0.7, and thus it is apparent that the slope 
of the linear rise is not the dominant factor determining 
the abundance of CO in average clouds. 

The plot with C4, the strength of the non- linear far-UV 
rise, (left panel of Figure 4 - same as bottom right panel of 
Figure 3) shows the CO abundance being correlated with 
stronger far-UV curvature. We note that for C4 > 0.5 we 
always detect CO and only upper limits exist for C4 be- 
low 0.25. This suggests that the carrier of the non-linear 
rise has a more direct association with the CO abundance 
than that of the linear rise; (c2)- Tlic right panel of Figure 
4 shows the relationship between the CO column density 
per E{B — V) and C4 assuming 6=1 km s~^ and T = 4 K. 
Although these values are likely not correct for each indi- 
vidual star it demonstrates the effective range of column 
densities observed and the lower limit of detectability by 
lUE. This correlation produces a Spearman's Rho of 0.630, 
consistent with the equivalent width relationship. 

The abundance of CO along a given sight line is gov- 
erned by the balance between formation through chemi- 
cal networks involving species such as C"*", OH, and CH, 
and destruction through photodissociation (van Dishoeck 
& Black 1988; Federman & Huntress 1989). The corre- 
lation of CO with C4 suggests that the presence of the 
carrier of the non-linear far-UV extinction is conducive to 
the presence of CO. This can be accomplished either by 
not destroying CO through shielding of the photodisso- 
ciating flux or by aiding in the formation of CO by not 
photodissociating the precursor molecules. 

CO photodissociates through the discrete absorptions of 
far-UV photons resulting in the excitation of the molecule 
into predissociating states. At the edges of clouds, this oc- 
curs primarily through the E — X (0-0) band at 1076 A. 
However, as optical depth increases into the cloud, most of 
the photodissociation moves to shorter wavelength transi- 
tions (van Dishoeck & Black 1988), where the extinction 
of flux by the non-linear rise increases dramatically, van 
Dishoeck & Black (1989) have demonstrated that the pho- 
todissociation rate of CO is reduced in cloud models that 
use steeper extinction curves to model the continuum ex- 
tinction of far-UV flux. It is therefore possible that we are 
observing an increase in molecular content associated with 
a decrease in photo destruction since the strength of the pa- 
rameters C2 and C4 regulate the far-UV flux at short wave- 
lengths. This is suggested in Figure 5, where we show the 
relation of these parameters to the CO normalized equiv- 
alent width. Howc;vcr. it is important to reahzc that the 
actual attenuation of radiation in a cloud depends strongly 
on the far-UV scattering characteristics of the grains re- 
sponsible for the extinction, even when the extinction rises 
sharply in the far-UV (Flannery et al. 1980). 



Although Figure 5 is suggestive, it is far from conclusive. 
If the carrier of the non-linear rise aids in the formation of 
CO, or any of its precursor molecules, such as CH and OH, 
the CO formation rate may be enhanced enough to explain 
the rapid increase in column densities seen in Figure 4. 



5. CONCLUSIONS 

The abundance of CO, as measured by the equivalent 
width of the UV absorption of the A — X (2-0) band nor- 
malized to miit E{B — V), is correlated with the strength 
of the non-linear far-UV rise in the UV extinction curves 
toward Galactic O and B stars. This correlation is in- 
dicative of either a decrease in the photo destruction rate 
of CO with increased extinction or of a dust environment 
that is conducive to the enhanced formation of CO or its 
precursor molecules. The linear rise does not correlate 
with CO abundance, although for sight lines with large Ry 
(> 3.6) we observe mostly small columns (< 10^^ cm^^ per 
E{B — V)). These results, together with those of Cardelli 
(1988) reinforces the notion that the use of the visual ex- 
tinction. Ay, as the independent variable in global studies 
of interstellar chemistry has to be viewed with some cau- 
tion. 

Future study of these relationships could benefit from a 
quantitative determination of the CO column density and 
an exploration of the H2/CO ratio in these environments. 
For H2, the ideal instrument is the recently launched Far- 
Ultraviolet Spectroscopic Explorer (FUSE) satellite. FUSE 
will be able to measure the H2 column from the individual 
lines of the Werner and Lyman band systems, as well as ob- 
serve the absorptions from the CO Hopfield-Birge bands. 
However, more accurate CO observations could be made 
by the Space Telescope Imaging Spectrograph (STIS) for 
two reasons. The wavelength separation of the rotational 
lines within the A — X bands is greater than that of the 
B — X and C — X and the STIS Echelle modes have higher 
resolution than FUSE, allowing for the measurement of the 
ground state rotational temperature and doppler velocity, 
and therefore a precise derivation of column density. Fi- 
nally, we note that the ratio of true absorption to scatter- 
ing by dust in the far-UV is critical to the interpretation of 
these results. Studies aimed at investigating this relation 
in conjunction with the FM extinction parameterization 
would be useful. Such measurements would facilitate a 
more detailed analysis of the conditions in these cloud en- 
vironments. 

All lUE data were obtained from the Astrophysics Data 
Facility at NASA's Goddard Space Fhght Center. The 
authors would like to thank Derek Massa for valuable dis- 
cussions. We also thank the referee, Steve Federman, for 
several helpful comments and suggestions. This work was 
supported by NASA grant NAG5-5122 to The Johns Hop- 
kins University. 



APPENDIX 

Typically, the column density is derived from the measured equivalent width through the use of an appropriate curve- 
of-growth. For molecules, the shape of a band integrated curve-of-growth depends strongly on the rotational excitation 
temperature, T, of the molecules and the line-broadening doppler parameter, h. Unfortunately, the spectral resolution 
of the lUE data is not high enough to allow for a direct measurement of the h parameter and only in the cases of the 
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strongest absorptions is it possible to put limits on the rotational temperature from the observed absorption profile. 

The use of multiple bands could produce a more accurate column density determination by limiting the possible values 
of b and T. However, to fit the relative amounts of absorption in multiple bands with a single model based on three 
parameters (A'^, b, and T) it is desirable to have accurate equivalent widths of more than three bands. The (0 0) and 
(1-0) bands of the A — X band system have blended -'^^CO and -'^''CO components. The (3-0) and (4-0) bands can 
sometimes be measured although the continuum has definite structure for some of the spectral types of stars considered in 
this study, while the (5-0) band is blended with the Si IV A1394 feature. Thus, the measurement of the equivalent width 
of these bands is dependent on an accurate fit to the continuum in these regions. The (6-0) band is not usually blended 
with other features, but it has an oscillator strength five times weaker than the (2-0) band and is only observed for the 
largest column density sight lines. Figure 6 shows band integrated curves-of-growth for the (2 0), (3-0), and (4-0) bands 
of the CO A — X band system computed using the band oscillator strengths and line wavelengths of Morton & Noreau 
(1994). The (4-0) band is blended with the a' — X (14-0) band due to quantum mechanical mixing of these states. All 
of these bands show the effects of saturation above the lUE detection limit. As can be appreciated from the following 
example, attempts to fit the equivalent widths of multiple bands can produce non-unique solutions to A'', b, and T. A 
typical measurement (HD 192281 e.g.) yields W(2-q) = 140 ± 13 ml, lF(..5_o) = 130 ± 16 mA, and T4^(4_o) = 94 ± 13 mA. 
Utilizing the curves-of-growth shown in Figure 6 we find that within the range of reasonable 6-values and rotational 
excitation temperatures the uncertainty in column density is approximately two orders of magnitude. 

Another approach is to synthesize a model spectrum that can be compared directly to the data. We attempted to solve 
for A'^, 6, and T by generating a grid of models with various values for these parameters and minimizing the value 
determined. This method proved successful only for those sight lines with the highest column density (~ 10"'^'' — 10^^ 
cm~^). The column density derived this way was usually within a factor of 2 or 3 of the valued determined by assuming 
T = 4 K and 6=1 km s~^. However, for smaller column sight lines (~ 10^^ cm~^) the column derived differed by as 
much as an order of magnitude. This is due to the fact that the absorption profiles for these columns at the resolution 
of lUE strongly mix the effects of differing T and b, especially at low T. These results suggest that the assumption of 
average values for T and b may be justified for a statistical comparison to the extinction curve parameters, but may not 
produce reliable answers on an individual basis. 
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Table 1 

Stellar extinction curve and CO absorption information 



Sp. Wco (2-0)^ 



Star 


Type 


E{B - V) 


C2 


C4 


(mA) 


Ref. 


SWP number(s)'' 


HD 13268 


08V 


0.44 


1.03 


0.38 


32 ± 23 


1 


09323 


HD 14434 


06V 


0.48 


0.80 


0.48 


< 67 


1 


16094 


HD 14442 


06IH 


0.73 


0.69 


0.44 


< 160 


1 


53141 


HD 15558 


05IH 


0.83 


0.59 


0.44 


106 ± 40 


1 


08322s 


HD 15629 


05V 


0.75 


0.63 


0.43 


135 ± 25 


1 


10754 


HD 21483 


B3III 


0.55 


0.730 


0.658 


316 ± 37 


2 


38298 


HD 34078 


09.5V 


0.53 


0.571 


0.520 


< 126 


2 


22108, 37429 


HD 36629 


B2V 


0.27 


0.21 


0.31 


< 42 


1 


10583 


HD 36879 


08IH 


0.50 


0.62 


0.28 


34 ± 13 


1 


02711s, 04654s, 04699s 


HD 36982 


B2V 


0.34 


0.050 


0.311 


< 36 


2 


13733, 21973 


HD 37022 


06p 


0.34 


0.033 


0.186 


< 16 


2 


07481 13737 13798 14597 15799 19606 


HD 37023 


B0.5V 


0.37 


-0.083 


0.153 


< 33 


2 


02770, 05017s 


HD 37061 


B0.5V 


0.54 


0.109 


0.044 


< 27 


2 


07038, 08086 


HD 37903 


B1.5V 


0.35 


0.384 


0.440 


< 23 


2 


06953, 08058, 21293, 29094 


HD 38087 


B5V 


0.33 


0.230 


0.311 


< 24 


2 


30046, 32748, 32749, 32750 


HD 46056 


08V 


0.51 


0.857 


0.541 


148 ± 37 


2 


08846s 


HD 46106 


BOV 


0.45 


0.591 


0.485 


55 ± 23 


2 


16435 


HD 46149 


08.5V 


0.48 


0.65 


0.59 


94 ± 28 


1 


21585, 28200 


HD 46150 


06V 


0.45 


0.65 


0.59 


127 ± 19 


1 


18947, 18948, 21584 


HD 46202 


09V 


0.47 


0.864 


0.515 


146 ± 32 


2 


30299 


HD 46223 


04V 


0.54 


0.70 


0.60 


107 ± 15 


1 


08138s 08338s 08844s 10757 


HD 47129 


08III 


0.36 


0.95 


0.46 


38 ± 10 


1 


10689, 13924 


HD 48099 


07V 


0.27 


0.874 


0.339 


< 20 


2 


52812, 52826, 52874 


HD 48279 


08V 


0.43 


0.72 


0.47 


60 ± 20 


1 


06504s 


HD 48434 


BOIII 


0.28 


1.10 


0.52 


< 46 


1 


06447s 


HD 73882 


08.5V 


0.72 


0.788 


0.540 


246 ± 18 


2 


30076, 30077, 38299, 38311 


HD 91824 


07V 


0.27 


0.633 


0.473 


< 87 


2 


16533 


HD 93028 


09V 


0.24 


0.811 


0.166 


< 69 


2 


05521 


HD 93204 


05V 


0.43 


0.48 


0.37 


< 32 


1 


07023, 07960 


HD 93205 


03V 


0.37 


0.45 


0.50 


< 28 


1 


09635s, 09655s, 09672s, 09738s 


HD 93222 


07III 


0.40 


0.626 


0.236 


< 76 


2 


22105 


HD 93250 


03V 


0.47 


0.58 


0.50 


70 ± 20 


1 


01618s, 06435s, 14746s 


HD 93403 


05III 


0.53 


0.61 


0.47 


92 ± 17 


1 


09075s, 09673s, 09739s 


HD 96715 


04V 


0.41 


0.67 


0.79 


198 ± 19 


1 


22000, 43980, 43981 


HD 147888 


B3V 


0.52 


0.133 


0.339 


38 ± 13 


2 


05159s, 05160s 


HD 147889 


B2V 


1.09 


0.151 


0.709 


221 ± 40 


2 


03924 


HD 147933 


B2V 


0.47 


0.139 


0.349 


79 ± 10 


2 


25428, 25429, 25430, 25431, 25439 


HD 149757 


09.5V 


0.32 


0.900 


0.563 


81 ± 8 


2 


36152, 36155, 36162, 36164, 36167, 38410, 41190 


HD 152233 


06III 


0.45 


0.74 


0.45 


93 ± 14 


1 


40981, 40987, 40988 


HD 154445 


BIV 


0.42 


0.309 


0.503 


< 44 


2 


06760s 


HD 162978 


08III 


0.35 


0.56 


0.17 


< 15 


1 


36945 36948 36949 38302 38303 


HD 164816 


BOV 


0.30 


0.37 


0.35 


< 25 


1 


02814s, 15308 


HD 165052 


07V 


0.43 


0.39 


0.28 


< 18 


1 


15306, 17106, 45466, 45502 


HD 167771 


07III 


0.44 


0.574 


0.453 


62 ± 13 


2 


09623s, 09633s 


HD 168076 


04V 


0.80 


0.48 


0.46 


127 ± 24 


1 


28277 


HD 192281 


05V 


0.70 


0.77 


0.43 


140 ± 13 


1 


43147, 44636, 44637, 46539 


HD 193322 


09V 


0.41 


0.879 


0.115 


< 15 


2 


31326, 31327, 31328, 31329, 36935, 36936 


HD 199579 


06V 


0.37 


0.898 


0.453 


139 ± 13 


2 


31317, 31318, 31323, 31331 


HD 200775 


B2V 


0.54 


0.54 


0.52 


276 ± 23 


1 


09836, 09837, 54089 


HD 215835 


06V 


0.65 


0.84 


0.57 


54 ± 34 


1 


45467, 45503 


HD 216532 


08V 


0.86 


0.53 


0.40 


171 ± 31 


1 


34226 


HD 216898 


08.5V 


0.85 


0.567 


0.333 


169 ± 34 


2 


43934 
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Table 1 — Continued 





Sp. 








Wco {2-0^ 






Star 


Typo 


E{B - V) 




C4 


(mA) 


Rcf. 


SWP nuniber(s)'^ 


HD 217979 


B2V 


0.61 


0.781 


0.580 


297 ± 92 


2 


43935 


HD 242908 


04V 


0.60 


0.54 


0.45 


< 58 


1 


16092 


HD 303067 


B0.5V 


0.31 


0.733 


0.147 


< 63 


2 


23759 


HD 303308 


03V 


0.45 


0.55 


0.43 


55 ± 22 


1 


09015s 


BD +60 594 


09V 


0.66 


0.819 


0.411 


103 ± 67 


2 


16644 


BD +60 2522 


BIV 


0.72 


0.66 


0.41 


< 114 


1 


08840 


CPD -59 2600 


06V 


0.53 


0.601 


0.179 


< 48 


2 


07021 



"errors are 1 a, upper limits are 2-a 

''an "s" indicates small aperture observation. 

References. — (1) Jenniskens & Greenberg 1993; (2) Fitzpatrick & Massa 1990 



Table 2 
asurv results 

Wx/E{B - V) 

versus Spearman p t value^ probability'' 





0.236 


1.83 


0.072 


7 


0.158 


1.21 


0.232 


Cl 


-0.369 


-3.00 


0.004 


C2 


0.395 


3.25 


0.002 


C3 


0.203 


1.56 


0.123 


C4 


0.667 


6.76 


8 X 10"' 



H = py jn^, where = 59 is the number of data points. 
''Probability that result is consistent with the null hypothesis. 
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Fig. 1. — Sample spectra showing an example of a strong CO A — X (2-0) absorption (HD 200775) and weak (HD 303308). The rest 
wavelengths of the R(0) Unes for the ^^CO and ^^CO are indicated. 
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Fig. 2. — Log-log plot of ^'^CO column density in cm ^ comparing our data to that listed in FedermaJi et al. (1994) for those sight lines 



in common. Column densities were determined using a curve-of-growth method assuming a doppler parameter b ■■ 
excitation temperature T = 4 K. Upper limits (2-(t) are shown as open triangles. 
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Fig. 3. — Correlation diagrams for the normalized equivalent width of the CO A — X (2-0) band in mA per magnitude of E{B — V). Upper 
limits (2-it) are shown as open triangles. 
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Fig. 4. — Left panel: Correlation diagram for the normalized equivalent width in mA per magnitude of E{B — V) versus C4. Upper limits 
(2-(t) are shown as open triangles. Right panel: Correlation diagram for column density in cm per magnitude assuming T = 4 K and b = 
1 km s~^. 
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Fig. 5. — Plot of C2 V. C4. The size of the circle indicates the strength of the CO absorption. Open circles are 2-a upper limits.The dotted 
lines correspond to constant k{x) for A = 1076 A, with increasing extinction to the upper right. 
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Fic;. 6.— Curvcs-of-growth for the CO A - X (2-0), (3-0), and (4-0) bands. The (4-0) band is blended with the a' - X (14-0) band. 
Tlic left ijaiicls are for 6 = 1 km s^-' with daslied, solid and dotted lines eorresponding to a rotational temperature of T = 10, 4, and 3 K 
respeetively. The right panels are for T = 4 K with dashed, solid and dotted lines corresponding to a doppler parameter b = 2, 1, and 0.5 km 
s^^ respectively. A horizontal dashed line shows our typical 30 mA detection limit. The top left panel shows a typical error bar in equivalent 
width. 



